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Background 



To investigate the molecular basis for invasion of malignant gliomas, proteomic analy- 
sis approach was carried out using two human glioma cell lines, U87MG and U343MG-A that demon- 
strate different motility and invasiveness in in vitro experiments. 



Methods 



High-resolution two-dimensional gel electrophoresis and matrix-assisted laser-de- 
sorption/ionization time-of-flight mass spectrometry analysis were performed. 

Nine distinct protein spots that were recognized with significant alteration between 
the two cell lines. Five of these protein spots were up-regulated in U87MG and four were up-regulated 
in U343MG-A. 



Conclusion 



Among these proteins, cathepsin D was shown to be one of the important proteins 
which are related with glioma invasion. However, further studies are necessary to reveal the exact role 
and mechanism of cathepsin D in glioma invasion. 
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INTRODUCTION 

Gliomas are a heterogenous group of neoplasms derived 
from various glial cells and account for 40-50% of all intra- 
cranial tumors. Glioblastoma multiforme (GBM), the most 
malignant type of glioma (WHO grade IV), is highly invasive 
and tend to diffusely infiltrate into the surrounding normal 
brain parenchyma. These characteristics prevent the tumor 
cells to reside within the margins of the therapeutic resection. 
Although the mechanisms of invasion in gliomas have been 
investigated in different human samples and experimental 
studies, there is a need for a novel panel of biomarkers that 
characterize each invasive phenotype. 

Traditionally, gene expression profiles have identified many 
genes with distinct expression patterns among different his- 
tological types and grades of gliomas [1-3]. However, biolog- 
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ical systems comprise protein components resulting from 
transcriptional and post-transcriptional control, post-transcrip- 
tional modifications, and shifts in proteins among the different 
cellular compartments. These properties cannot be analyzed 
by microarray systems at the RNA level, whereas proteomic 
analysis such as high-resolution two-dimensional gel electro- 
phoresis (2DGE) and mass spectrometry (MS) separate and 
identify the targeted proteins [4,5]. Proteomics has recently 
emerged as a new field of protein science that allows for high- 
throughput profiling of expressed proteins at cellular and sub- 
cellular levels. This proteomic -based analysis promises new in- 
sights into the cellular mechanisms involved in tumor invasion 
and is likely to result in the discovery of novel diagnostic and 
therapeutic targets. 

In our previous study, different migration ability between 
human GBM cell lines was confirmed by simple scratch mi- 
gration assay [6]. Furthermore, we identified some genes 
possibly related to the high motility in GBM cell lines using 
differentially display-polymerase chain reaction (PCR) by re- 
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verse transcription-PCR and Northern blot analysis [6], and 
have continued to reveal the downstream mechanism for these 
genes [7-10]. 

In this study, based on the migration/invasion ability of the 
high motile cell line (U87MG) and less motile cell line (U343- 
MG-A), we examined these cell lines using a proteomics ap- 
proach, including the use of 2DGE and matrix assisted laser- 
desorption/ionization time-of-flight mass spectrometry 
(MALDI-TOFMS) analysis, to compare the proteomic profiles 
of these two cell lines with different motility 

MATERIALS AND METHODS 

Cell lines and culture conditions 

U87MG and U343MG-A (the American Type Culture 
Collection, Rockville, MS, USA) was maintained in Dulbec- 
co's modified Eagle's medium (DMEM, Gibco, Carlsbad, CA, 
USA) media supplemented with 10% fetal bovine serum (Gib- 
co, Carlsbad, CA, USA) under 5% C0 2 at 37°C 

Matrigel invasion assay 

2X104 cells of U87MG and U343MG-A in 0.35 mL of se- 
rum-free DMEM were seeded into the upper chamber of a 
10- well chemotaxic chamber (Neuro Probe, Gaifherburg, MD, 
USA) and complete DMEM was placed in the lower cham- 
ber, and a Matrigel-coated membrane was inserted between 
the two chambers. After 24 hours incubation at 37°C, mem- 
branes were fixed and stained using a hemacolor rapid stain- 
ing kit (Merck, Darmstadt, Germany), the cells from 5 ran- 
dom microscopic fields (220x magnification) were counted. 
Same experiment was repeated 3 times, independently. 

Scratch test 

For complete inhibition of cell proliferation, 5 mM hydroxy- 
urea was added to the culture media for 24 hours. The cul- 
tures were then scratched with a single-edged razor blade, 
washed twice with PBS, and placed in media containing hy- 
droxyurea. After 46 hours incubation, the cultures were im- 
mersed in methanol and stained with 0.1% toluidine blue. 
Three microscopic fields were evaluated for each wound inju- 
ry. The number of cells that traversed the scratched area and 
the maximum distance traveled were determined in each field 
and averaged for each injury. These experiment were repeated 
three times. 

Preparation of proteins from cell lines 

For proteome analysis, both cell lines were passed when an 
approximate 80% confluence was reached, then each cell line 
was homogenized in lysis buffer consisting of 7 M urea, 2 M 
thiourea, 2% CHAPS, 1% difhiofhreitol (DTT), 0.5 mM EDTA, 



and 1 x protease inhibitors. The samples were centrifuged at 
14000 g for 20 min at 18-20°C Supernatants were then sepa- 
rated and used for 2DE after protein quantification using the 
BSA assay (Bio-Rad Laboratories, Munich, Germany). 

Two-dimensional gel electrophoresis 

The first dimension of 2DE was performed on an immobi- 
lized pH gradient (IPG) strip (Amersham Biosciences, Amer- 
sham PI, UK). Linear pH 3~10 IPG strips (18 cm) were rehy- 
drated overnight at room temperature in rehydrating buffer (8 
M urea, 1% DTT, 2% CHAPS, and 0.5% IPG buffer). Sample 
of 300 mg was applied during rehydration. The first dimen- 
sion was run for 53500 Vh at 20°C using the following condi- 
tions: 500 V for 1 h, 1000 V for 1 h; and 8000 V for 6 h and 30 
min. Next, gels were equilibrated for 30 min in a equilibration 
buffer I containing 50 mM Tris-Cl, 6 M urea, 30% glycerol, 2% 
sodium dodesyl sulfate (SDS) and 0.1% DTT and subsequent- 
ly equilibration buffer II containing 50 mM Tris-Cl, 6 M urea, 
30% glycerol, 2% SDS and 0.25% iodoacetic acid. The second 
dimension was run on Ettan DALT II system (Amersham Bio- 
sciences, Amersham PI, UK). An 8-15% gradient acrylamide 
gel was used for the second dimensional gel electrophoresis. 
The IPG strips were placed on the surface of the second di- 
mensional gel, and then were sealed with 0.5% agarose in SDS 
electrophoresis buffer containing 25 mM Tris base, 192 mM 
glycine and 0.1% SDS. The gels were placed into Ettan DALT 
II system chamber containing 1% SDS electrophoresis buffer. 
The gels were run overnight at 100 V until the dye front reached 
the bottom of the gel. 

Proteins were visualized using coomassie blue staining, in 
which the intensity of the protein spot is correlated with the 
quantity of protein present. The gels were fixed at room tem- 
perature in 12% tri-chloroacetic acid for 60 min and washed 
twice in distilled water for 15 min. Next, the gels were staining 
overnight with 160 mL of staining solution (0.1% coomassie 
G-250 in 2% H 3 P0 4 , 10% ammonium sulfate) plus 40 mL of 
methanol, and washed for 1 to 3 min in 0.1 mol/L Tris, H3PO4 
buffer, ph 6.5, then rinse for 1 min in 25% aqueous methanol. 
The gels were then stabilized the protein-dye complex in 20% 
aqueous ammonium sulfate. 

Image analysis and sample preparation 
for mass spectrometry 

Coomassie blue-stained gels were analyzed to select spots 
that appeared consistently between two samples. Selection of 
the spots was based on volume and density as determined 
using Phoretix 2D Advanced software (Version 5.01; Nonlin- 
ear Dynamic, Newcastle, UK). The selected proteins were ex- 
cised under a laminar flow hood using a scalpel and then pro- 
cessed for protein digestion. An accelerated digestion protocol 
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was adopted as described by Havlis et aL [11]. Briefly, the Coo- 
massie blue-stained protein spots were excised, destained, and 
dehydrated in ammonium bicarbonate and acetonitrile (ACN) 
for 20 min. The solution was then aspirated, and gel pieces were 
brought to complete dryness in a SpeedVac. Gel pieces were re- 
hydrated with a freshly prepared solution of sequence-grade 
modified trypsin (20 ng/uL in 25 mM ammonium bicarbonate 
buffer) at 4°C for 60 min and then incubated at 55°C for 60 min. 
The tryptic peptides were extracted with 10 uL of 1% formic acid 
aqueous solution by sonication in a water bath for 15 min. The 
resulting samples were stored at -80°C until analyzed by MS. 



peptides and sequence coverage: five or more matched peptides 
or 30% or greater sequence coverage was considered acceptable. 

Statistical analysis 

Quantitative data were collected for invasion assay, scratch 
test, and migration assay, as described above. These experi- 
ments were performed in triplicate and repeated three times. 
The observed effects was evaluated by the Mann-Whitney U 
test using SPSS version 20.0 software program for Windows 
(SPSS Inc., Chicago, IL, USA). The level of significance was 
set atp<0.05. 



MALDI-TOFMS and protein identification 

Tryptic-digest peptides were concentrated and desalted us- 
ing GELoader microcolumns packed with POROS 20 R2 
resin (Applied Biosystems, Foster City, CA, USA). The pep- 
tides were then eluted with 1 (±L matrix solution (70% ACN, 
0.1% trifluoroacetic acid). The eluted peptides were then 
spotted on to a stainless steel target plate and peptide masses 
were determined using MALDI-TOFMS and an Applied 
Biosystems 4700 proteomics analyzer (Applied Biosystems). 
The unit was calibrated using external mass standards 
(Calmix 1 and 2; Applied Biosystems). Peptide masses were 
compared with the theoretical peptide masses of all human 
proteins in the NCBI database using Mascot software (Matrix 
Sciences, London, UK). A Mascot probability-based score for 
protein identification of greater than 50 was considered accept- 
able. The identified proteins were also examined for matched 



RESULTS 

The different motility between U87MG 
andU343MG-A 

In invasion assay, the number of invasive cells (mean± 
standard deviation) was 179.3±15.1 in U87MG and 25.5±7.5 
in U343MG-A. There was a significant difference in both cell 
lines (p<0.001), as shown in Fig. 1. The number of migrated 
cells was 473±112.4 in U87MG and 80.2±13.5 in U343MG- 
A. This difference was also statistically significant (p<0.001), 
as shown in Fig. 2. 

Protein expression patterns and differential analysis 
of U87MG and U343MG-A 

Protein expression maps of U87MG and U343MG-A were 
generated using 2DGE. The 2DGE gel images of these cell lines 
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Fig. 1. Matrigel invasion assay. The number of invasive cells (mean±SD) was 179.3±15.1 in U87MG (A) and 25.5±7.5 in U343MG-A (B). 
U87MG cell line was more invasive than U343MG-Awith a statistical significance (p<0.001) (C). SD: standard deviation. 
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Fig. 2. Scratch test. The number of migrated cells (mean±SD) was 473±112.4 in U87MG (A) and 80.2±13.5 in U343MG-A (B). U87MG cell 
line was more motile than U343MG-Awith a statistical significance (p<0.001) (C). SD: standard deviation. 
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demonstrated similar patterns of protein expression, due to the 
genetic similarity between them (Fig. 3). With Phoretix 2D Ad- 
vanced software for 2DGE gel analysis, numerous well- stained 
protein spots were detected. Only nine spots revealed significant 
difference in spot density and volume. Compared to the spot in- 
tensity and volume of U343MG-A, spots 509, 561, 566, 613, 618 
were up-regulated and spots 171, 274, 342, 475 were down-reg- 
ulated in U87MG cells. The result of different expressed protein 
spots was summarized with circles and spot numbers as shown 
in Fig. 4. 



Identification of protein spots using 
MALDI-TOFMS 

We analyzed differentially expressed protein spots, followed 
by excision from the gel, using MALDI-TOFMS and database 
searching. The 9 identified proteins in spots are summarized in 
Table 1. The higher score or greater number of matched pep- 
tide/sequence coverage indicated the better probability. Pro- 
teins that were up-regulated in U87MG included one protein- 
ase protein (cafhepsin D), two structural proteins (eukaryotic 
translation inhibition factor 6 and nicotinamide phosphoribo- 
syl transferase), one RNA-binding protein (ribonuclease in- 
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Fig. 3. Coomassie blue-stained two-dimensional gel electrophoresis gel images of U343MG-A and U87MG. 
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Fig. 4. Protein spots differentially expressed between U343MG-A and U87MG. Note that the spot intensity & volume of U87MG are described 
compared to those of U343MG-A. 
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hibitor 1) and one metabolic enzyme protein (transglutamin- 
ase 2). Down-regulated proteins in U87MG were one metabolic 
enzyme protein (glyceraldehyde-3-phosphate dehydrogenase), 
one ATP-binding protein (autophagy related 7), one F-box pro- 
tein (F-box only protein 47) and one signal transduction pro- 
tein (calcium/calmodulin-dependent phosphodiesterase 1A). 

DISCUSSION 

The current study describes the proteome maps of human 
GBM cell lines, U87MG and U343MG-A. Considering the fact 
that U87MG is a more motile cell line than U343MG-A, some 
of the different expressed proteins may be related with motil- 
ity or invasion in glioma cell lines. Two cell lines showed simi- 
larities in the protein expression patterns on the 2DGE gel im- 
ages, reflecting that there is closeness in the genetic background. 
This point was observed in another proteomic analysis using 
sibling glioma cell lines, J3T-1 and J3T-2, with different inva- 
sion property [12]. Among the only 18 protein spots that are 
differently expressed and identified, annexin A2 was focused 
on as a candidate invasion-related protein, involved in angio- 
genesis-dependent patterns in gliomas. 

In our study, one of the proteins identified, cathepsin D, 
was expressed at higher levels in U87MG than in U343MG-A. 
Cathepsin D is normally localized within the lysosomes and 
involved in degradation of proteins and processing of precur- 
sor proteins [13]. It plays an essential role in tumor progres- 
sion, angiogenesis and apoptosis, in addition to the protein 
turnover process [14-16]. In a recent investigation, cathepsin 
D prevented oxidative stress-induced cell death through the 
activation of autophagy in cancer cells. This result suggests that 
inhibition of autophagy may be a novel therapeutic strategy 
for human cancers which express cathepsin D [17]. Fukuda 
et al. [18] suggested that cathepsin D may be a potential se- 
rum marker for the prediction of aggressive nature of human 
gliomas. In neuroblastomas, cathepsin D has an important role 
in anti-apoptotic signaling and is a novel mechanism for the 
chemo-resistance [15]. 

Based on proteomic analysis, Iwadate et al. [19] found that 
37 proteins differentially expressed based on the histologic 
grading of human gliomas. Among these proteins, cathepsin 
D was increased in high-grade gliomas and may be one of 
novel biomarkers for survival prediction and rational target for 
anti-glioma therapy. Recent investigations revealed that ca- 
thepsin D was confirmed to be up-regulated in high-grade 
classification and possibly related to acquired temozolomide 
resistance in human malignant gliomas [20]. Liu et al. [21] pro- 
posed that the inhibition of lysosome cathepsin D exocytosis 
from glioma cells plays a pivotal modulatory role in their mi- 
gration and invasion. Furthermore, promotion of prolifera- 



tion and invasion, and suppression of apoptosis in glioma cells 
induced by Rab27a were exerted by increasing cathepsin D and 
regulated by miR-24 [22]. Our study also implicates the possi- 
ble role of cathepsin D in glioma motility and invasion. How- 
ever, further studies using knock-down or stable transfection 
are necessary to reveal the exact role and mechanism of cathep- 
sin D in glioma invasion. 

Transglutaminase 2, one of the up-regulated proteins in 
U87MG in our experiment, has not been widely investigated 
in glioma research. However, recent papers proposed that 
this protein is possibly related to drug resistance in glioma 
cells [23]. Although there have been some papers showing 
the inhibitory role of ribonuclease inhibitor 1 on cancer growth 
and metastasis [24], there have been no past studies address- 
ing the role of ribonuclease inhibitor 1 in gliomas. Nicotin- 
amide phosphoribosyltransferase has been known as one of 
the key enzymes in salvaging pathway for the synthesis of 
nicotinamide adenine dinucleotide, involved in cell metabo- 
lism and proliferation [25]. Considering the high metabolic 
demand in malignant gliomas, the inhibition of nicotinamide 
phosphoribosyltransferase may reduce nicotinamide adenine 
dinucleotide to inhibit the growth of glioblastomas [26]. 
However, there is no previous study on the relation between 
nicotinamide phosphoribosyltransferase and glioma invasion. 
This has not been investigated in the role of eukaryotic trans- 
lation inhibition factor 6 in gliomas. 

In conclusion, using proteomic analysis for two glioma cell 
lines with different motility and invasiveness, some proteins 
were differently expressed. Although there was no protein di- 
rectly connected to tumor cell motility among them, cathep- 
sin D seemed be one of the possible proteins which are relat- 
ed to glioma invasion. 
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